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We study structure formation in cosmological void regions using high-resolution 
hydrodynamical simulations. Despite being significantly underdense, voids are popu- 
lated abundantly with small dark matter halos which should appear as dwarf galaxies if 
' their star formation is not suppressed significantly. We here investigate to which extent 

, the cosmological UV-background reduces the baryon content of dwarf galaxies, and 

thereby limits their cooling and star formation rates. Assuming a Haardt & Madau UV- 
' background with reionisation at redshift z — our samples of simulated galaxies show 

i that halos with masses below a characteristic mass of Mf.{z = 0) = 6.5 x 10^ h^^ Mq 

' are baryon-poor, but in general not completely empty, because baryons that are in 

fH \ the condensed cold phase or are already locked up in stars resist evaporation. In halos 

Q^. with mass M < Mc, we find that photo-heating suppresses further cooling of gas. The 

I ' redshift and UV-background dependent characteristic mass Mc{z) can be understood 

O , from the equilibrium temperature between heating and cooling at a characteristic 

" ■ overdensity of (5 ~ 1000. If a halo is massive enough to compress gas to this density 

despite the presence of UV-background radiation, gas is free to 'enter' the condensed 
phase and cooling continues in the halo, otherwise it stalls. By analysing the mass 
accretion histories of dwarf galaxies in voids, we show that they can build up a sig- 
nificant amount of condensed mass at early times before the epoch of reionisation. 
Later on, the amount of mass in this phase remains roughly constant, but the masses 
' of the dark matter halos continue to increase. Consequently, photo-heating leads to a 

reduced baryon fraction in void dwarf galaxies, endows them with a rather old stel- 
lar population, but still allows late star formation to some extent. We estimate the 
resulting stellar mass function for void galaxies. While the number of galaxies at the 
faint end is significantly reduced due to photo-heating, additional physical feedback 
processes may be required to explain the apparent paucity of dwarfs in observations 
of voids. 
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1 INTRODUCTION 

Large regions of space that contain few or no galaxies can 
be clearly identified in modern spectroscopic redshift sur- 
veys. Ab out 25 years ago, such 'voids ' have been first dis- 
covered llGregB QL fc Thompsonl Il978l: iJoeveer et J] llQTSt 
iKirshner et al.lll98 j) . but it remains a challenge to explain 
why they are apparently so empty. 

It is well known that hierarchical models of struc- 
ture formation in standard cold dark matter (CDM) cos- 
mologies produce large underdense regions in the distri- 
bution of matter jPeebledllQsS iHoffman fc Shaham|[l983 



Ivan de Weveaert fc van KampenI Hool) . but these under- 
dense regions still contain structural elements and bound 
halos, even though the characteristic masses of these ob- 
jects are several orders of magnitude smaller than corre- 
sponding ones found in average regions of the universe. 
fCottlobcr ct a!L (200^) predicted that a typical 20 Mpc 
diameter void should contain up to 1000 halos with mass 
~ 10^ M© and still about 50 halos with mass ~ 
10^" /i"^ Mq. Assuming a magnitude of Ms = —16.5 for 
the galaxy hosted by a halo of mass 3.6 X lO"/i-^M0 
liMathis fc Whit&,2002f) predict that about five such gala^c- 
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ies should be found in the inner regions of a typical void of 
diameter 20 Mpc. 

Over the last decad e, there were many attempt s to find 
dwar f galaxies in voids Lin^i£^^^^96; Podcs cu et alJ 
ll997l:lKuhn et al.lll997l: iGroein fc Gelleil »1999). An overall 
conclusion from these studies has been that faint galax- 
ies do n ot tend to fill u p the voids outlined by the bright 
galaxies. IPeeblej ll200lll pointed out that the dwarf galax- 
ies in the Optical Redshift Survey (ORS) follow the dis- 
tribution of bright g alaxies remarkably clo sely. Using the 
SDSS data release 2, iGoldberg et al.l (|200^ measured the 
mass function of galaxies that reside in underdense regions. 
They selected galaxies as void members if they had less than 
three neighbours in a sphere with radius 7 Mpc. More 
than one thousand gala:xies passed their selection criterion 
(which differs slightly from the criterion used in the numeri- 
cal simulations of Gottlobcr ct al. (2003) ^-^d this paper). 
Their measurements are consistent with the predictions from 
the numerical simulations once the tendency of more mas- 
sive halos to concentrate at the outer parts of voids (where 
they still may pass the nearest neighbour selection criterion) 
is taken into account. However, the observational situation 
is unclear for halo masses smaller than ~ 10^''/i~^Mq. At 
present, there are no observational hints that a huge num- 
ber of dwarf galaxies in voids may exist, despite the large 
number of small halos predicted by CDM models. 

This findin g resembles the 'substructure problem ' in 
galactic halos (jKlvpin et al.l Il99fll : iMoore et alJ Il999^ . A 
solution for both problems could arise from physical pro- 
cesses capable of suppressing star formation in dwarf galax- 
ies. The two major effects proposed in this context are su- 
pernova feedback and heating of the gas in halos by the 
UV-background radiation. The latter increases the ther- 
mal pressure, and as a result, the gas in systems with 
10^ K can be evaporat ed out of halos 



memura fc Ikeuchilll984l : iDekel fc ReesHlQSTl : lEfstathioul 
ri99a: iBabul fc ReesI Il992l) . Similarly, supernova feedback 



could drive a galactic outflow that remo ves a signifi- 
cant fraction of the gas in a dwarf galaxy jPekel fc Silk! 
[l986; C ouchman fc Rees. .l986). However, the efficiency of 
supernova driven winds depends strongly on the details 
of model assumptions |Navarro fc Steinmet j Il997l . I2OOOI: 
ISprineel fc Hernauistll2003^ . A general result is that low- 
mass systems are much more e asily affected by supernova 
driven winds than larger ones iMac Low fc Ferraral [19991 
iNavarro fc Steinmetzll200ol) . 

During the epoch of reionisation, the gas temperature 
is raised to a few times 10* K. ^ees ( 1986) argued that in 
dark matter halos with virial velocities around ~ 30 kms~^ 
gas can then be confined in a stable fashion, neither able 
to escape nor able to settle to the centre by cooling. The 
gas in significantly smaller systems is thought to b e virtu- 
ally e vaporated [Thoul fc Woinbere 1996; Barkana fc Loebl 
Il999tl . However, due to central self-shielding the evapora- 
tion o f gas that has already coo led can be significantly de- 
layed JSusa fc Umemurall2004al[bl). Us ing radiative transfer 
simulations. ISusE^^^n^muraT 2004lJl pointed out that ha- 
los may contain a significant amount of stellar mass pro- 
duced before reionisation occurred, even if their remaining 
gas mass is evaporated during reionisation. 

In a 3D Euleri an adaptive mesh refinement simulation, 
iTassis et alJ ll2003l) found that the global star formation 



rate was significantly reduced during the reionisation epoch. 
Their result also indicated that stellar feedback enhances 
this effect dramatically. Thus, an imprint of the epoch of 
reionisation may be expected for the stellar population in 
dwarf galaxies. Indeed, almost all dwarf galaxies app ear to 
have an early epoch of star formation (lMatedll998l) . but 
there is no distinct t ime at which star formati on becomes 
generally suppressed jGrebel fc Gallaghe J l2004fl . 

In this paper, we use high-resolution hydrodynamical 
simulations of cosmological void regions to analyse the star 
formation and cooling processes of void galaxies. In partic- 
ular, the simulations are well suited for studying the evolu- 
tion of isolated dwarf galaxies from the epoch of reionisation 
to the present. This allows us to examine whether the UV 
background has a sufficiently strong effect on dwarf galax- 
ies, keeping them faint enough such that their abundance 
can be reconciled with observations. We also determine the 
characteristic mass scale below which cooling is suppressed 
by the UV background. An analysis of the spatial distribu- 
tion of dwarfs, the impact of supernova feedback, and the 
spectral properties of the stellar content of the formed dwarf 
galaxies will be discussed separately. 

Our study is organised as follows. Details of our simu- 
lations are described in Section |5| In Section |3| we analyse 
first the baryon fraction as a function of both halo mass and 
redshift. Then we investigate the mass growth of the con- 
densed phase which consists of both stars and cold dense 
baryons. We identify the characteristic mass below which 
halos are subject to evaporation. Finally, we estimate the 
galaxy mass function in void regions. We discuss and sum- 
marise our results in Sectional 



2 SIMULATIONS 
2.1 Numerical method 

Our simulations have been run with an updated version o f 
the parallel Tree-SPH code Gadget JSpringel et alJl200lll . 
The code uses an entropy- conserving formulation of SPH 
JSprineel fc Hernauistll2002^ which alleviates problems due 
to numerical over-cooling. The code also employs a new algo- 
rithm based on the Tree-PM method for the A'^-body calcu- 
lations which speeds up the gravitational force computation 
significantly compared with a pure tree algorithm. 

Radiative cooling processes for an optically thin pri- 
mordial mix of helium and hydrogen are included, as well 
as photo-ionisation by an external, spatially uniform UV- 
background. The net change of thermal energ y content is 
calcul ated by closely following the procedure of iKatz et alJ 
il996l) for solving the rate equations. We also use the set of 
cross sections cited in this paper, but adopt a slightly mod- 
ified Haardt fc Madau ( 1996) UV-background with heating 
rates as depicted in Fig.Q 

The physics of star formation is treated in the code by 
means of a sub-resolution model in which the gas of the inter- 
stellar medium (ISM) is described as a multiphase medium 
of hot and cold gas (Yep es ct al . 1992,; . Spi;ingeLfc_Hcr na uis3 
2003) . Cold gas clouds are generated due to cooling and are 
the material out of which stars can be formed in regions 
that are sufficiently dense. Supernova feedback heats the 
hot phase of the ISM and evaporates cold clouds, thereby 
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Table 1. Main characteristics of the void simulations. A/gas and A/jark denote t he mass of a gas and o f a dar k matter particle in the 
simulation, respectively. Feedback parameters are according to model described in lSpringel fc Hernguisj j20o3). The UV-flux, Jq^ , at 
2 = is given in units of 10~^^ ergss"^ cm~^ sr~^ Hz~^. For the description of the simulation with an additional heat pulse see Sec. l3.2( 
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Figure 1. Photo-heating rates due to the ambient UV- 
background as a function of redshift. We adopt a time evolu- 
tion for the UV-background as given bv lHaardt fc Madaul ri996ll . 
Reionisation in this model takes place at 2 = 6. 

establishing a self-regulation cycle for star formation. The 
heat input due the supernovae also leads to a net pressuri- 
sation of the ISM, such that its eftective equation of state 
becomes stiffer than isothermal, see Fig. |21 This stabilises 
the dense star forming gas in galaxies against further gravi- 
tational collapse, and allows converged numerical results for 
star formation even at moderate resolution. We also follow 
chemical enrichment associated with star formation, but we 
have neglecte d metal-line cooling in compu ting the cooling 
function. See ISprineel fc Hernauisd l)2nn3|) for a more de- 
tailed description of the star formation model implemented 
in the Gadget code. 



2.2 Initial conditions 

Using the mass refinement technique described by 
iKlvpin et alJ (|200J^ . we simulate void regions with high 
mass resolution, embedded in a proper cosmological en- 



vironment. Our voids have been selected for resimulation 
from two periodic computational boxes of side-lengths L — 
80h~^ Mpc and 50h^^ Mpc, respectively. To construct suit- 
able initial conditions, we first created an unconstrained ran- 
dom realization at very high resolution, using the ACDM 
power spectrum of perturbations. For the large box, A'^ — 
1024'^ particles were used, while for the smaller box, we em- 
ployed 2048^ (~ 8.6 billion) particles. The initial displace- 
ments and velocities of the particles were calculated using 
all waves ranging from the fundamental mode k = 2n/L to 
the Nyquist frequency kny = 2n/L x N^^^ /2. To produce 
initial conditions at lower resolution than this basic high- 
resolution particle set-up, we then merged particles, assign- 
ing to merged particles a velocity and a displacement equal 
to the average values of the original small-mass particles. 

In this way, we first run 128^ low-resolution simulations 
until the present epoch and selected the void regions from 
them. T he algorithm for identify ing the voids is described in 
detail in lGottlober et al] i2003r) . It allows us to select void 
regions with arbitrary shape. To this end, the method starts 
from a spherical representation of the void which is then ex- 
tended by spheres of smaller radius, which are added from 
the surface of the void into all possible directions. However, 
in the present application we have restricted the resimula- 
tion to a spherical void region to avoid ambiguities in the 
definition of allowed deviations from spherical shape. 

In the second simulation step we use the original sam- 
ple of small-mass particles in the regions of interest when 
we construct initial conditions. Thus we reach a mass reso- 
lution within the void regions that corresponds to the 1024'' 
or 2048^ set -up, respectively. We use a series of shells around 
the voids where we progressively merge more and more of 
the particles until the efi'ective resolution of 128^ particles 
is reached again far away from the voids. This procedure 
ensures that the voids evolve in the proper cosmological en- 
vironment and with the right gravitational tidal fields. Mix- 
ing of particles of different mass occurs only in the shells 
surrounding the high resolution voids. Finally, we split the 
particles in the regions of high mass resolution into dark 
matter and gas particles. For all simulations, we adopted a 
concordance cosmological model with = 0.3, Q,a = 0.7, 
= 0.04, h = _ffo/(100kms-i Mpc-^) = 0.7 and as = 0.9. 
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Figure 2. The distribution of particles in the density - temperature phase diagram for different redshifts for the high-res run. The 
soUd curves give the equihbrium hne. It is derived from the heating-cooling module in the code by computing the temperature at which 
the heating rate equals the net cooling rate. We assume that the 'entry' into the condensed branch occurs at 10^ (Pbaryon)- The circles 
indicate the derived entry temperatures for the individual redshifts. See Sec. lij.3i and Fig. llUl for a discussion of the equilibrium and entry 
temperatures. 



2.3 Simulation runs 

Using the multi-mass technique described above, we have 
re-simulated a void region in the 80 Mpc box with three 
levels of refinement. The mass of a dark matter particle in 
the void is 3.4 x IO^/i^^Mq. The corresponding SPH gas 
particles have an initial mass of 5.5 X IO'^^-^Mq. Note that 
some gas particles may reduce their mass during the run 
(or vanish entirely) if they undergo star formation and cre- 
ate new coUisionless star particles. In our analysis of the 
simulation results, we in general only consider halos com- 
posed of a minimum of 150 dark matter particles. For the 
80 Mpc simulations this corresponds to a minimum halo 
mass of 5.1 x 10® h^^ M© and a circular velocity ~ 23 km/s. 

We have also carried out re- simulations of a void re- 
gion in the 50 Mpc box. This leads to a substantially 
improved mass resolution with the same level of refinement. 
Here we achieve a mass resolution of 8.2 x 10^ /i"^ Mq for 
the dark matter particles, corresponding to a minium halo 
mass of 1.2 X 10® /i"^ M0. In addition, we have evolved this 
region also with the full resolution available based on the ini- 
tial high-resolution particle set-up (the 2048^ particle grid, 
corresponding to four levels of refinement) . In this case, the 
meiss resolution is improved to 1.0 x 10® M© for the dark 
matter particles, and the minium halo mass reaches down to 
1.6 X lO*/i-iM0. We give an overview of our simulations in 
Table [J where w e also list the main simulation parameters. 

iPower et alj i2003l) gave a simple criterion for the grav- 
itational softening length necessary in A^-body simulations, 
e <; ''200 / V N200 ■ Obeying this condition ensures that par- 
ticles do not suffer stochastic scattering in the periphery 
of t he halo which exceed the mean gravitational accelera- 
tion. IPower et alJ ll2003ll found that even a larger softening 
length keeps the central density profile unaffected. However, 
in a hydrodynamical simulation including radiative cooling 
we wish to use a softening as small as possible to obtain an 
optimal resolution for the cooled gas. Our halos are in the 
mass range from about 10® Mq to 10^^ Mq. Using 

4 1 

M200 = 200 X -Trr^oo (p), 
where (p) is the average cosmic matter density, we find the 



corresponding softening lengths, e ~ 2.2 to 1 kpc. Since 
we carry out simulations including radiative cooling we lower 
the softening slightly and use for all simulations the maxi- 
mum between 2 kpc comoving and 0.8 kpc physical. 
For the very high-resolution run, the parameters lh~^ kpc 
and 0.5 h^^ kpc, respectively, are used. We have imposed a 
minimum SPH smoothing length equal to the gravitational 
softening length. 



In order to analyse in more detail the effects of the 
UV-background on the baryonic content of halos, we have 
also carried out several additional runs of our basic simu- 
lation of the void identified in the 50 h"^ Mpc box. Here 
we used different choices for the star formation and feed- 
back parameters, and three different values of the UV-flux 
normalisation, spanning four orders of magnitude. We also 
ran the same simulation without thermal stellar feedback. 
Using this run we can demonstrate that thermal feedback 
itself has only a minor impact on the halo baryon fraction. 
For all simulations with star formation and feedback, we se- 
lected similar parameters fo r the mul tiphase model of the 
ISM as used bv lSpringel fc HernauistI (.2003.) . However, we 
here have not included kinetic feedback (wind model) from 
supernova, since then it would be difficult to disentangle the 
effects caused by the UV background from those caused by 
supernova driven winds. 



The simulations were performed on parallel supercom- 
puters, an IBM Regatta p690 (Jiilich Supercomputer Cen- 
ter, Germany), a SGI Ahix 3700 (CIEMAT, Spain) and on 
several different Beowulf PC clusters at the AIP and lU 
Bremen. The typical CPU time for a simulation with up to 
5 million particles was ~ 9 CPU days on an SGI ALTIX 
with 32 processors. The highest resolution simulation with 
44 million particles was run on an AMD Opteron Beowulf 
cluster and took a little less than 2 months of CPU time 
using 64 processors. 
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Figure 3. Mass function n(> M) for our basic run halo sample 
at 2 = 0. The mass function is derived by taking into account 
the total mass within the virial radius. The small vertical line 
indicates a mass of 150 dark matter particles, which we consider 
as a lower limit for an acceptable resolution. The numbers along 
the line indicate the actual number of halos. The dashed line is 
obtained by considering the dark matter mass instead of the total 
mass. For comparison the modifie d Sheth-Tormen mass function 
derived bv lGottlober et al.l j20o3) is shown (dotted line). For our 
halo sample, the circular velocity Hcirc = GMtot/rvir can be well 
approximated by Vdrc = 31 kms~i X (Mtot/lO^O h'^ Mq)°-^'^. 



3 SIMULATED DWARF GALAXIES IN VOIDS 
3.1 Halo mass function 

We identify virialised halos usinR th e Bound Density Max- 
ima (EDM) algorithm jKlvpin et al.l lToQg'l. In this method, 
galaxy halos are found from local density maxima with an 
iterative procedure to identify the centre of mass in a small 
sphere around the centre. Then, radial density profiles are 
computed. Particles that are not gravitationally bound to 
the system are excluded in the computation of the total 
mass. The radius of the system is selected as the minimum 
between the virial radius and the point at which the density 
profile stops declining (e.g. because a nearby halo is encoun- 
tered, or the halo lies within another halo). We define the 
virial radius as the radius where the enclosed mean density 
equals the value expected for a top-hat collapse model, 



4/37r) 



(1) 



where (p) is the mean cosmic matter density. For the case of 
a flat cosmology with Qm + = 1, a useful approximation 
for the red shift dependent characte ristic virial overdensity 
IS given bv lBrvan fc NormanI Jl998l) : 
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Figure 4. Baryon fraction in individual halos for differently re- 
solved simulations. We compute for each halo the baryon fraction 
within the virial radius, (M* -|- Mgas)/(Mt + Mgas -I- A^dm)- We 
take only those halos into account which consist of more than 150 
dark matter particles. From the left to the right the arrows indi- 
cate the smallest resolved halos in the high-res, basic, and void2 
run, respectively. For the high-res run we approximate the baryon 
fraction by Eq. J3 (solid line). 
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with the cosmological expansion factor a = 1/(2; -I- 1). The 
mean density (p) evolves with redshift as 



{p){z) = p< 
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In low density regions, the number of interacting halos or ha- 
los with substructure is very small. Thus, the radii of virtu- 
ally all our objects correspond to their spherical-overdensity 
virial radii. For the same reason, the fraction of unbound 
particles in the halos is small. Hence, we can simply consider 
all particles within the virial radius to compute further halo 
properties. 

In Figure El we show the cumulative mass function 
for our sample of void halos, based on the measured to- 
tal virial masses. For comparison, we also include a line 
for the dark matter halo mass function alone, which how- 
ever show s appreciable differenc es only for the most mas- 
sive halos. ICottlober et al.l (JioO^) derived a modified Sheth- 
Tormen mass function for halo populations in underdense 
regions. They showed that the mass function can be derived 
from the mean density, flm,void, in the volume considered. 
For our halo sample analysed here (based on the high-res 
run), the relevant mean density amounts to Jlm.void — 0.03. 
The predicted mass function for this value is in good agree- 
ment with the one measured for the simulated halo sample. 
Our simulations thus are in good agreement with conclusion 
obtained in previous works: Voids are filled with a significant 
number of halos with masses M < 10^° /i-i Mq. If each of 
these halos contains the mean cosmic baryon fraction, and 
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cooling and star formation within them is not suppressed sig- 
nificantly, a high density of luminous dwarf galaxies should 
be expected in voids. 

3.2 Baryon content of dwarf galaxies 

In Figure 3] we show the baryonic mass fraction, = 
Mbaryon(< Tvir ) /Mtot ( < ^vir), for each halo identified in 
our simulated voids. The more massive halos in our sam- 
ple, Aftot 2 X 10^" Mq, have approximately the cos- 
mic mean baryon fraction, fh,cosm = ^b/^m- However, for 
smaller halos, the baryon fraction decreases rapidly with de- 
creasing halo mass. In fact, most of the smallest halos are 
nearly free of baryons. 

We find that the baryon fraction is insensitive to the 
details of the definition of the virial radius of halos, because 
the cumulative baryon fraction varies only very slowly in 
the outskirts of halos. This can be seen in Figure |K| where 
we show the radial profile of the baryon fraction for halos 
of different mass. While the sizes of halos can be systemati- 
cally affected by different definitions of the virial radius (one 
may for example choose to use only the dark matter for the 
definition and not the total mass), the measured baryonic 
fractions are robust. 

We quantify the transition between the two extremes, 
'baryon-rich' and 'b a ryon- poor', using the fitting function 
proposed bv lGnediiil (l200(jl . 

h = hn U + {-2''"-l)i^Xy"\ (5) 



where /^o is the asymptotic baryon fraction in massive halos. 
Figure |S| indicates that a baryon fraction decreasing with 
mass can be reasonably well approximated assuming a = 2 
in Eq. O. lGnedinl i2000fl found a less steep transition, which 
may be caused by radiative transfer effects included in his 
code. A partial self-shielding in halos may then reduce the 
radiative heating. For very small halos the approximation 
seems to fail in any case. We discuss a possible origin for 
this in Sec. 13.51 

In the very massive halos of our samples we found a 
roughly constant baryon fraction fbo, independent of red- 
shift and numerical resolution. The value of fbo ~ 0.16 we 
measured lies slightly above the cosmic mean. At the char- 
acteristic mass, Mc, the baryon fraction is /b = /6o/2 by 
definition. For 2 = 0, we derive a characteristic mass of 
Mc = 6.5 X 10^ Mq from the 'high-res' run. 

It is important to consider whether numerical resolution 
effects influence this result. In particular, a too small num- 
ber of resolution elements per halo could easily introduce a 
spurious baryon reduction due to numerical oversmoothing. 
The resulting characteristic mass would then depend on the 
number of particles in a halo rather than on the halo mass 
itself when simulations with quite different mass resolutions 
are compared. In Fig. 2] we compare the void2, basic and 
high-res runs. Their baryon fractions as a function of halo 
mass overlap nicely, even though the mass resolution differs 
by more than an order of magnitude. Therefore, significant 
numerical oversmoothing occurs only in halos smaller than 
those included in our analysis. 

The marginal displacement of the characteristic mass 
as a function of resolution can be understood in terms of 
the better resolved merging histories of higher resolution 
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Figure 5. The baryon fraction as a function of radius for halos 
with different mass. All halos are chosen at 2 = from the basic 
simulation. Radii are normalised to the BDM virial radius. Labels 
along the lines indicate the mass of the halos. 



runs. The better the mass resolution, the more small pro- 
genitors of a final halo can be resolved and contribute stars 
to the final object. In particular, many more star particles 
are formed while the UV-background radiation is still van- 
ishingly small. The resulting difference can be seen in Tab. 1: 
The basic and the high-res run have identical initial condi- 
tions, except for the increased resolution, which leads to the 
production of ~ 10% more stars. As a result, the characteris- 
tic mass of the high-res run is also slightly lowered. However, 
this effect is limited by the fraction of mass in progenitors 
which can form stars before reionisation. 

The importance of the UV-background for the baryon 
fraction is demonstrated in Fig.Q We have carried out three 
simulations for which we multiplied our standard UV-flux 
for every redshift by the factors 0.0, 0.01, and 100, keep- 
ing the spectral shape of the UV-background radiation con- 
stant. These strong variations of the UV-background flux 
level displace the characteristic mass-scale only by a factor 
< 3. We will discuss in the next section. Sec 13.31 that this 
shift is caused by the modification of the equilibrium tem- 
perature. The run with zero UV-flux clearly shows that the 
background radiation causes the baryon deficit: In this case, 
even the smallest halos show the average cosmic baryon frac- 
tion. 

During the epoch of reionisation, the true heat input 
may be larger than computed in our heating-cooling scheme, 
which is based on the assumption of coUisional ionisation 
equilibrium. In the onset of reionisation, non-equilibrium ef- 
fects can be significant, however, and the ionised fraction 
increases very rapidly. While our implicit solver for the evo- 
lution of the thermal energy deals gracefully with this situa- 
tion, part of the injected energy may be missed if the ionized 
fraction jumps from zero to a finite value in the course of 
one time step. One may speculate that the heat pulse at 
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Figure 6. The baryon mass as a function of the total mass en- 
closed in the virial radius from the high-res run. Solid squares 
indicate halos with stars, and crosses those without. The dotted 
line shows the dependence if each halo would contain the mean 
cosmic baryon fraction. The solid line is derived based on the 
assumption that in low-mass halos the gas density follows the 
average distribution in the p-T diagram (see Sec. l3.5i for more de- 
tails). The arrow indicates a baryon mass consisting of 150 SPH 
particles. The properties of smaller halos may be affected by the 
poor SPH resolution in these halos. Thin and thick dash-dotted 
lines show the approximation of Eq. j^J, with exponents a = 1 
and 2, respectively. 

reionisation is underestimated by our method, but that it 
may have a significant effect on the baryon fraction in the 
hafos fater on if fully taken into account. In order to test this 
question, we have carried out two simulations which mimic 
an upper limit for the impulsive heat input: from z — 7 to 
6, all gas is heated to a minium temperature of 1.5 x 10 K 
and to 5.0 x 10* K in the runs imp-heat-15 and imp-heat- 
50, respectively. After z=6 the gas evolves again according 
to our standard heating-cooling scheme, see Fig. |H| Even 
with a strong heat input at the epoch of reionisation, our 
measured characteristic mass at z = is hardly aflected, 
as can be seen in Fig.[7| Moreover, the different heat input 
at reionisation has virtually no effect on the evolution of 
the characteristic mass. This is caused by the short cooling 
times at high redshift; Fully ionised gas with average cosmic 
density and a temperature of 5.0 x 10'' K has at 2: = 6 a 
cooling time as short as ~ 50Myr, because the maximum of 
hydrogen line cooling is in this temperature range. Inverse 
Compton cooling with CMB photons is also very efficient at 
these redshifts. It keeps the cooling times short even if we 
would heat to a temperature above the efffcient line cool- 
ing. Since gas in the halo of a galaxy is much denser, the 
cooling time is even shorter. Hence, the energy injected by a 
heat-pulse at the epoch of reionisation is radiated away on 
time-scales smaller than the dynamical ones. 

Finally, efficient stellar feedback can in principle also 
remove gas from small halos. However, for the thermal feed- 
back considered in the multiphase feedback model used in 




Figure 7. Baryon fraction as a function of mass for different UV- 
fluxes. The photo-heating by the UV-background is varied from 
zero to = 95 X 10"-^'' ergs s"^ cm~^ sr"^ Hz"^ , see Tab. 1. 

The average baryon fraction in mass bins is computed; for the 
basic-run also the r.m.s. deviation is shown by error bars. In ad- 
dition baryon fractions are depicted for the impulsive heat model 
with 1.5 X 10* K (imp heat) and for the simulation without ther- 
mal stellar feedback (no-feed). 



our simulations, such a gas removal does not occur. While 
the feedback regulates the consumption of cold gas by star 
formation, it does not cause gaseous outffows. The latter 
only occur in our simulations if explicitly modeled with a 
kinetic feedback component (Springol & Hornauist 200^. 
However, we deliberately avoided the inclusion of such feed- 
back models in this study, allowing us to focus on the impact 
of UV- heating in a clean fashion. Figure [7| verifies that in 
simulations without feedback the characteristic mass scale 
is at the same place as in the basic run. 

In order to determine the evolution of characteristic 
mass with redshift, we use a least-square fit of Eq. @ to 
our measurements of the baryon mass fraction at a num- 
ber of different simulation output times. We can infer the 
transition mass scale reasonably well from our high-res sim- 
ulation up to 2 ~ 5. In Fig. 1^1 we show the resulting evolu- 
tion of Afc(z). Interestingly, as structure grows towards lower 
redshifts, progressively more massive halos become baryon- 
poor. 

It is interesting to also compare this result with the 
time evolution of Ale measured for the simulation with an 
additional heat pulse at the epoch of reionisation, which is 
also shown in Figure IHl The two curves are identical ex- 
cept for a small offset. As discussed above, the latter can be 
attributed to the different mass resolutions of the two sim- 
ulations. The similarity of the two curves therefore shows 
that the heat pulse at reionisation does at most weakly in- 
ffuence the characteristic mass at 2: < 3. Mc{z) is relatively 
insensitive at these redshifts to the previous thermal history 
of the low density gas at high redshift, but it depends more 
strongly on the current state of the gas and that of the UV 
background. 
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Figure 8. The distribution of particles in the density - temperature phase diagram for different redshifts for the impulse-heat-15 
simulation. From z = 7 to 6 the minimum gas temperature is set to 1.5 X 10'' K. For smaller redshift the standard heating-cooling scheme 
is applied. At z ~ 4 the initial heat pulse is faded away by adiabatic cooling for densities Pgas/(p}baryon ^ O-l- Equilibrium and entry 
temperatures are indicated by the solid lines and open circles, respectively. See Sec. I3.3l and Fig. llUl for a discussion of them. 



The evolution of the characteristic mass Mc{z) can be 
expressed as 



Mc(z) 



IQio /i-i Mp 



r{z) 



1 



1 



3/2 



Ac(0) 



1/2 



(6) 



were t{z) encodes the evolution of the minimum virial tem- 
perature required for halos to still allow further cooling in 
the presence of the UV-background. We will discuss this cri- 
terion and its derivation in more detail below. Treating t(z) 
as a simple analytic fitting function for the moment, we find 
that our numerical results can be well described by 



t{z) = 0.73 x{z + if-^^ exp{-(0.25 zf'^} 



(7) 



for the 'high-res' run over the entire redshift range, see 
Fig.H 



3.3 How to suppress gas condensation 

A quantitative understanding of the characteristic mass 
can be obtained by considering the equilibrium line be- 
tween photo-heating rate and cooling rate in the density- 
temperature phase-space plane. In Figure \^ we show the 
gas particles at three different times, and include this equi- 
librium line for high-density gas, computed self-consistently 
from the cooling and heating routines in Gadget. Much of 



10^ 



10 is indeed dis- 



the gas in the density range 5 
tributed along this line. 

We define the equilibrium temperature at a fiducial 
overdensity of 1000 as "entry temperature" Tcntry into the 
condensed phase. Gas that condenses in a halo will at least 
reach this temperature due to photo-heating before it can 
cool further to the ~ 10* K reached at very high overden- 
sities. Note that this is independent of the potential differ- 
ence between a "c old mode" or a "hot mode" of accretion 
jKeres et al J 120051) . In the cold mode, gas creeps along the 
lowest possible temperature into the condensed phase, with- 
out being heated by an accretion shock to the virial temper- 
ature first, as it happens in the "hot mode". However, even 
in the cold mode, the gas will at least be heated to the "entry 
temperature" Tontry by the UV-background. After reaching 
this temperature, it can then evolve along the equilibrium 



line towards higher densities and eventually reach the onset 
of star formation. In a sense, for halos with Tvir = Tcntry, the 
hot and cold mode should therefore become largely identical. 

We argue that a comparison of Tvir with Tontry provides 
a simple criterion that tells us whether the gas in a halo can 
still cool. To demonstrate this, we show that this assumption 
provides a quantitative explanation for our measurements of 
Mc{z). We define the virial temperature for our halos as 



fee Tvir 



1 GMvir 
-fJ-TUp 



(8) 



where mp is the proton mass and fca is the Boltzmann con- 
stant. The mean molecular weight of the fully ionised gas 
is /i = 0.59. Virial mass and radius are related by the def- 
initions of Eqs.Q and (|1J. Hence the virial temperature 
depends on the mass as 



_ 1 finip 



Ac(2)n„ 



{l + z){GMv.rHo}i 



(9) 



For the cosmological concordance model used here we insert 
f^m = 0.3 and cast Eq. (|5J into 



Ac(0) 



M ■ 



T.r = 3.5xlO*K(l + .)^^^^^^^ I 1010,- 

(10) 

Now we apply our criterion introduced above: If the virial 
temperature of a halo is below the entry temperature into 
the condensed phase, its potential well is not deep enough 
to compress the gas sufficiently and to overcome the pres- 
sure barrier generated by the photo-heating. Thus, we ex- 
pect that the characteristic halo mass necessary to increase 
the condensed baryonic mass can be estimated by rewriting 
Eq. (Cnj into 



Mc(z) 



1010 /i-Wf 



T 



entry 



1 



3.5 X 104 K 1 + z 



A,(0) 
A,(z) 



(11) 



This relation can only be an approximation, since Mc{z) will 
in general depend on the entire history of a halo, whereas 
Tontry is au iustantaueous quantity. However, we have seen 
from the 'imp-heat' simulation that the characteristic mass 
is mainly determined by the current state of the gas. We 
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Figure 9. Evolution of the cliaracteristic mass scale Mc for sev- 
eral simulations. We approximate the halo baryon fractions for 
each simulation by Eq. J3 ■ We determine Mc using a least-squares 
method. ff,Q is a free parameter, but it is almost constant for 
all simulations, /i,o ~ 0.16. Open circles indicate the results for 
the high-res simulations. Here the smallest halos included have 
1.5 X 10* Mq. Crosses and squares show the results for imp- 
heat-15 and imp-heat-50 simulations, respectively. Here the small- 
est halos included have 10^ Mq. Note that the characteristic 
mass may lie somewhat below the mass of the smallest halos, 
since it is derived from a fit to the baryon fractions. The solid 
line shows the approximation for Mc{z) from the high-res run 
using Eq. J^J and Q. 



therefore assume that Mc{z) can be considered as an in- 
stantaneous quantity. In this sense, Eq. IIIH gives the moti- 
vation for our ansatz of Eq. ^ . Also it shows that in a first 
approximation the fitting function t{z) is proportional to 



t entry 



In Figures |5| and |H| solid lines indicate the equilibrium 
temperature, Tl(Tcq, p) — r(req, p), for 5 > 10^. In that den- 
sity region cooling times are short and the gas is essentially 
distributed along the equilibrium line -or according to the 
multi-phase ISM prescription used in Gadget-2. The equi- 
librium temperature depends on the UV-background. The 
dependency on the amplitude of the background radiation, 
Jo, is shown in Fig. IIUI If we lower or increase the UV-flux 
by two orders of magnitude, the entry temperature shifts 
roughly by a factor of two. With Mc oc Tc^^ry ^6 obtain 
the shift of Ale as found in the simulations, see Sec. 13. 21 and 
Fig.lZI 

In the following, we use the abbreviation Tontry ~ 
TcntTy{z) /3.5 X 10* K. In Figure fTTI we show the evolution 
of rentry(z) for the high-res run and compare it with the 
expression for t{z) derived from fitting the measured red- 
shift evolution of the characteristic mass in the same sim- 
ulation. First of all, we note that the redshift variation of 
both is small compared with the evolution of the charac- 
teristic mass. This indicates that the large variation of the 
latter is mainly governed by the 1/(2-1-1) term in Eq. Q. 
Furthermore, from Fig. Illl we can conclude that the simple 
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Figure 10. The equilibrium temperature, Teq, as a function of 
the amplitude, Jq, of the UV-background radiation. Tcq is com- 
puted by solving 7i{Teq,p) = r{Teq, p). Heating and cooling func- 
tions as implemented for the numerical simulations are used. The 
amplitude dependence is given for three densities, p = <5f2bPcriti 
for z = 0. 



criterion invoked above reproduces the evolution of Mc{z) 
up to 2 ~ 5 astonishingly well, with a deviation of a few 
10%. 

However, the condition Tvir = Tentry appears to under- 
estimate the characteristic mass for small redshift to some 
extent. This is not surprising for a number of reasons. For 
example, we compute the virial radius from the dark mat- 
ter distribution only, hence the total masses used here are 
slightly above the virial masses. Furthermore, Tcntry(z) de- 
termined at a fixed overdensity 1000 may underestimate the 
true entry temperature. We can account for this and improve 
our match of the numerical results by refining our criterion 
in the following way: The virial temperature of a halo has 



to be > 1.3 X Tontry to permit further condensation. 

To explain the deviation at high redshifts we have to ac- 
knowledge the fact discussed above: t{z) and Tcntry are com- 
pletely different in the sense that the latter denotes which 
halos are able to add more gas to the condensed phase at 
a given time, while the former reflects the entire evolution 
history of dwarf halos. Before reionisation, no halo is photo- 
evaporated, i.e. Mc = r = 0. It takes some time after reioni- 
sation to produce baryon-poor halos. Hence, t{z) should be 



below Tt 



entry 



for some time after reionisation. 



In any case, the agreement demonstrated in Figure [TTI 
shows that the argument that the virial temperature should 
at least equal the entry temperature provides a good quan- 
titative description for the transition scale between baryon- 
poor and baryon-rich halos. Interestingly, this explanation 
can account for the effect simply by alluding to the ongoing 
photo-heating of the gas by the UV-background radiation, 
and the accompanying increase of the gas pressure in low 
mass halos. At 2 < 3 the effects of the impulsive heating 
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Figure 11. Evolution of the entry temperature. We derive from 
the high-res run the entry temperature, Tcntry, for several red- 
shifts (open triangles). Tentry is determined according to the de- 
scription in the caption of Fig. |^ For comparison we plot the 
expression t{z) given in Eq. Q (solid line). The scaled entry 
temperature ("refined model") is also shown (dashed line). 



during the epoch of reionisation play only a subdominant 
role. 



3.4 Baryon deficit in the linear theory 

Small-scale baryonic fluctuations grow slower than the cor- 
responding dark matter fluctuations due to the counteract- 
ing pr essure gradients. To describe this effect . ICnedin fc Hui 
l|l998r introduced a flltering wavenu mber fcp over which 
baryonic fluctuations are smoothed out. lGnedinI ll200Cf) com- 
pared the corresponding filtering mass 



My 



4tv , , { 2-Ka\ 



(12) 



up to 2: = 4 with the characteristic mass measured in cosmo- 
logi cal reionisation simula tions and found a good agreement. 

iGnedin fc Huil il998l) showed that fep can be related to 
the linear growth function D{t) by 



dt 



, cl {D + 2HD) 
47rG{p) 



dt" 



(13) 



with the time-dependent Hubble parameter H — d/a. We 
rewrite Eq. I13II in order to integrate it up to z = 0. The 
linear growth function obeys the relation 



D + 2HD = 4tvG{p) D, 



(14) 



which we can use to remove the time derivatives in Eq. 11311 . 
We also substitute the integration variable by a. With 



. Ho , 



= S\a) = l + !^„ 



i-l) +!^A(a'-f 1) (15) 



we then get 
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(16) 



For a given c osmology the linea r growth function can be 
computed by llCarroll et alJll992l) 



D{a) = -nr. 



S{a) 



da 



SHa) 



(17) 



hence we obtain a solution for Eq. 1131 for any redshift, 
provided the evolution of the sound speed, Cs(a), is known. 
To obtain the sound speed for t he high-resolution regi on in 
the simulations we follow again iQ nedin fc Hull Jl998l) and 
compute the volume averaged temperature. 



s(a) 



5 fee (r) vol (a) 
3 /imp 



(18) 



Figure 1121 shows that {T)voi rises sharply at the time of 
reionisation, thereafter it decreases slowly from ~ 5000 K 
to 1000 K. These temperatures are lower than expected for 
an average cosmic volume because the considered void re- 
gion is underdense and the temperature scales with density 
(see Fig.OJ. 

In Fig. 1121 we show the resulting filtering mass scale. For 
redshift z = 0, it is about three times larger than the charac- 
teristic mass obtained from the simulations, while the differ- 
ence becomes progressively smaller towards higher redshift. 
One may argue that the underdensity of the void region has 
to be taken into account. To consider a rather extreme case 
we assume that the void evolves like an open universe with 
f^m = 0.03 and unchanged f^A = 0.7. As shown in Fig. 1121 
(dashed line) the density of the region has only a moder- 
ate effect on the filtering mass. It reduces to some extent 
the difference with our measured characteristic masses. At a 
first glance one may expect a different density dependence of 
Mp. The average density enters Eq. 1121 directly via its ap- 
pearance in the denominator of Eq. 11311 . If the density did 
not affect the Hubble constant, H, and the linear growth 
function, D, we would expect the filtering mass is propor- 
tional to (p)^^^^. However, H and D depend on (p) as can 
be seen from Eq. 1141 . In order to derive consistently the 
filtering mass as a function of flm the density dependence 
of H and D, as given in Eqs. 1151 and 1171 . has to be taken 
into account. For high redshift, a ^ 1, we can approximate 

1 



This leads to 



and hence 



Mp oc n; 



-1/2 



For high redshift we obtain the same proportionality as ex- 
pected without taking the contribution of H and D into ac- 
count. However, the result changes for lower redshifts. For 
a = 1 we have S = 1, i.e. there is no dependence on flm- 
This implies that the integral Eq. 1161 depends less strong 
on f2,„ and in Mp the factor (p) becomes dominant. 

In summary, we expect that at high redshift the filter- 



ing mass shows the proportionality Mp oc ft, 



-1/2 



whereas 



for low redshift even the opposite behavior may occur. The 
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Figure 12. Comparison of filtering and characteristic mass. The 
thin fine indicates the evolution of the volume averaged temper- 
ature in the simulation. The thick solid line gives the filtering 
mass integrated according to Eq. I16i using the averaged tem- 
perature depicted here. The dashed line shows the filtering mass 
for a low density region with Qm = 0.03. Open squares indicate 
the characteristic mass obtained from the high-res simulation. 



actual dependence can only be found by computing the in- 
tegral in Eq. 1131 . Note that the speed of sound Cs is also 
a function of redshift and modifies the integral additionally. 
The dashed line in Fig. 1121 gives the result for a low den- 
sity region with Qm = 0.03 and for the rather low average 
temperature in the void region. At low redshift the filter- 
ing mass overpredicts significantly the characteristic mass 
we measure. 

The density dependence in Eq. H12|l is partially compen- 
sated by the fact that for a <^ 1 the filtering wave number 
scales with density according to fc| oc Qm- In conclusion, we 
find that at low redshift the filtering mass Mf overpredicts 
the characteristic masses we measure. 



3.5 Are there dark halos free of baryons? 

We argued above that photo-heating may prevent further 
gas condensation in dwarf halos. If it never takes place in a 
given halo, or an existing condensed phase has been evapo- 
rated at some time by the UV-background, the halo will only 
contain diffuse gas, where "diffuse" for our purposes denotes 
gas of too low density to support star formation. It is there- 
fore also interesting to distinguish between halos with and 
without stars, as we have done in Fig. |S] This allows us to 
demonstrate again how closely the equilibrium line in the p- 
T diagram, the virial temperature, and the baryon fraction 
are connected. 

The expected baryonic mass in a halo is given by 

Mb = ^vrrJi.Pb, (19) 

where p^, denotes the average gas density in a halo. At red- 
shift z = 0, the majority of the gas is distributed at low 



densities along the power-law relation 

[j^J ^ 3.6 X 103 K (2°) 

in the p-T phase diagram, where photo-heating balances adi- 
abatic cooling due to the expansion of the universe. If we 
assume that the average temperature in halos is given by 
the virial temperature, then the average density (of the dif- 
fuse gas) in a halo cannot lie below the relation given by 
Eq. II2UII . i.e. this imposes a rough upper limit on the av- 
erage gas density, and hence the diffuse baryonic content of 
the halo. With the help of Eq. IIOII . we can replace the virial 
temperature by the virial mass and cast this condition for 
the average gas density in a halo into 

(21) 

Plugging this result into Eq. 11911 and using the definition of 
the virial radius, Eq. iQ, we obtain 

Mb _ 84.7 Qb f Mvi, 1"'^ 
IQio/i-iMo Ac(0) nm llOio/i-iMoJ ■ ^ ' 

for the baryon mass at jz = 0. The resulting baryon mass 
is shown in Fig. |H] Halos without any stellar matter are 
primarily distributed just above this line. Hence Eq. 1221 
provides a rough upper limit for the diffuse baryonic mass 
in dwarf-sized halos. In summary, for halos which contain 
virtually no condensed baryons, the diffuse gas mass can be 
determined from the equilibrium line in the p-T phase dia- 
gram. This can hence explain why halos significantly below 
Ale are seemingly almost empty of gas. It provides also a 
lower limit for the baryon fraction. 

3.6 Condensation history 

The basic argument above is that photo-ionisation can stop 
at some time the condensation process of gas in dwarf galaxy 
halos. We can support this idea further by finding halos with 
a constant amount of mass in the condensed phase. Fieure lTHl 
shows the evolution of the different mass components for two 
example halos. A halo selected with mass significantly above 
Mc at z = has continuously increased the dark matter, 
stellar, and condensed masses, while the gas mass decreased. 
In a halo with total mass significantly below Mc, the indi- 
vidual mass components change as well, but in contrast to 
the case of the more massive halo, the condensed mass is 
remarkably constant. This lends strong support to our basic 
argument. We see here that for the individual halo shown 
neither further gas condensates nor photo-heating evapo- 
rates all the gas. Thus even in the optically thin treatment 
of the galactic gas reservoir, photo-heating is not able to ex- 
pel gas from the galaxy that has overdensities above ~ 1000. 
Instead, the condensed gas continues to be slowly converted 
into stars. It is thus not expected that photo-heating instan- 
taneously switches off star formation. 

We now analyse the mass accretion histories of several 
halos of different mass. In Fig. 1141 we show that the evo- 
lution of the condensed phase changes systematically from 
halo masses below Mc to those above. For the latter, the to- 
tal mass and the condensed mass grow almost in parallel. For 
halos with mass close to Mc, the condensed mass increases 
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Figure 13. Evolution of the mass in different components for two high-resolution halos, one above Mc (left panel) and one below Mc 
(right panel). The evolution of the total mass (thick solid line), the condensed gas mass (thick dashed), the stellar mass (thick dotted) 
and diffuse gas mass (thin dotted) is shown. In addition the evolution of the condensed mass (gas+stars) is shown (thin solid line). 



monotonically but slower than the total mass. Finally, for 
halos with mass significantly below Mc, the condensed mass 
remains constant. Thus the evolution of individual halos is 
perfectly consistent with the result that photo-heating pri- 
marily stops condensation in small halos. Moreover, we find 
that our "refined" criterion for condensation reproduces the 
time at which condensation stops quite well. To show this 
explicitly we have plotted in the left panel of Fig. 1141 Mc 
derived form Tvir = 1.3 x Tentry The times at which the 
low-mass halos fall below Mc correspond to the times at 
which the condensed masses begin to remain constant (right 
panel) . 

In halos with Mtot(O) < 2 x IO^'/i'^Mq the condensed 
gas phase is apparently evaporated completely. Even if these 
halo appear to be sufficiently well resolved at z = 0, we cau- 
tion that condensation can only take place at 2: < 4, prob- 
ably at even higher redshift since M[z) < Mc{z). Hence we 
expect that in a simulation with higher resolution the con- 
densed phase might stay and result in a noticeably larger 
amount of stars. However, we cannot really exclude that 
condensed gas in halos with Mtot(O) < 2 x IO^H-'^Mq be- 
comes generally unstable at some time and is therefore evap- 
orated even for arbitrary good numerical resolution. Such an 
assumption is frequently invoked i n simple analytic treat- 
ments. For instance, the analysis of lBarkana fc Loebl (Il999l) 
is based on the assumption that all gas with T > Tvir is evap- 
orated. The fact that the significantly improved resolution 
of the 'high-res' run compared to the 'basic' run does not 
lead to a stable condensed gas phase in all halos with mass 
about Mtot(O) ~ 2 x W^h~^MQ seems to point into this 
direction. However, more stringent test of this will require 
better resolution for the progenitor of these halos. 



3.7 Stellar mass function 

Using our estimates for the suppression of baryon condensa- 
tion in halos with mass below Mc(2:), an approximation for 
the expected stellar content in small halos may be derived. 
To this end we assume a mean mass accretion history for 
the small halo s in ou r sample. IWechsler et al.l i2002l) and 
Ivan den BoschI (|20o3) showed that the average dark matter 
mass accretion histories can be approximated by 

MMAH(z) = Mtot(2 = 0)e-"", (23) 

which works especially well for halos with an early forma- 
tion time. In Fig. 1141 we show an average accretion history 
(thin solid line ) using a = 0.35. Moreover, from the work of 
van den BoschI ll2002l) one can derive an appropriate value 
for a for halos with ~ 5 x 10^ Mq in a low density 
universe. 

All halos in our sample have rather flat accretion histo- 
ries compared to the evolution of the characteristic mass 
Mc{z), thus most of them have been able to condensate 
baryons and to form stars at high redshift. But at some time 
the halo mass fell below Mc{z), and from this time onwards 
the mass in the condensed phase is expected to remain con- 
stant. To test this, we estimate the stellar mass in a halo by 
computing first the redshift z®'' at which the mass accord- 
ing to the accretion history equals the characteristic mass, 
Mma-ii{z'"^) = Mc{z'"^). At this redshift, we determine the 
condensed mass, assuming here that at this time the baryon 
fraction in the halos has the cosmic mean value Qi,/Qrn- In 
addition, we assume that on average 80 % of the baryons 
are in the condensed phase, see Fig. 1131 Finally, as an upper 
limit for the amount of stars that may form, we assume that 
all of the condensed baryons are eventually converted into 
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Figure 14. Mass accretion histories. We compute the mass accretion histories by choosing from the high-res run several halos at z = 
and searching repeatedly for the most massive progenitor. The left panel shows the evolution of the total mass, while the right panel 
shows the evolution of the condensed baryon mass of the same halos. Jumps in the condensed mass indicate merger events. The numbers 
indicate the dark matter mass at 2 = in 10^" Mq. In the left panel, the evolution of the characteristic mass for the 'high-res' run is 
also shown (thick solid line). For comparison, we also plot the mean mass accretion history (thin solid line) for Aftot = 1.4 X 10^ Mq 
with a = 0.35, as expected for a low density environment, see Eq. ^231 . 



stars, (/cq) ~ 0.8. Hence we estimate the stellar mass in a 
given halo as Moq(2; = 0) = (/eq) Qb/^m Mmah(z°''). For 
halos more massive than Mc{z = 0) the maximum stellar 
mass is given just by (/cq) fib/fim Aftot(O). 



We define the stellar mass function, 7Vhaio(Af*), as the 
number of halos with stellar mass larger than M*, inde- 
pendent of the dark matter masses of the halos that host 
the galaxies. Figure IT^ shows the stellar mass function ob- 
tained by the procedure described above (solid line). For 
comparison, we also show the mass function for a fixed stel- 
lar meiss fraction /, for all halos (dashed line). The de- 
viation of the two curves below M, ~ 10^ Mq indi- 
cates that in low-mass halos the stellar mass fraction falls 
below the fixed value /, , hence those halos are shifted to 
smaller M* and the mass function begins to fiatten. For a 
stellar mass content of ~ 10* Mq the number of halos 
is reduced by a factor of three. The simulated halo sam- 
ple (dash-dotted line) provides a lower limit for the stellar 
mass function because of the limited numerical resolution. 
Above ~ 10* Mq, the two curves agree very well. For 
smaller masses, the two curves deviate since some halos have 
formed too few stars due to insufficient resolution. In con- 
clusion, the photo-evaporation reduces the number of halos 
with Mr > 10* Mq by a factor of about four and, more 
importantly, the faint-end slope of the stellar mass function 
is much shallower than expected based on the dark halo 
mass function of A'^-body simulations alone. 



4 DISCUSSION 

The main purpose of this work has been to investigate the 
effect of UV-heating on the baryon fraction and the stellar 
mass of dwarf-sized isolated halos in voids. This is moti- 
vated by the suggestion that photo-heating may provide a 
potential explanation for the observed low number density 
of galaxies in voids, and the apparent paucity of luminous 
galactic satellites in the Local Group, in contrast to the large 
abundance of dark matter halos and subhalos, respectively, 
predicted by CDM models. 

In order to examine this question, we have carried out 
some of the presently best resolved hydrodynamical simula- 
tions of galaxy formation in large regions of space. Our sim- 
ulations are unique in the sense that they specifically sample 
cosmological void regions at high resolution and evolve them 
to the present epoch. 

We have identified a characteristic mass scale below 
which the baryon fraction in a halo is reduced by photo- 
heating. This happens for masses M^iriz = 0) < 6.5 x 
10^ Mq (denoted here as dwarf galaxies), and this char- 
acteristic mass scale depends only weakly on the UV-fiux. 
The suppression itself happens due to photo-evaporation of 
gas out of halos, thereby reducing the baryon fraction of 
dwarf galaxies, and also by offsetting cooling losses in ha- 
los such that further condensation of baryons is strongly 
reduced. We have derived a simple quantitative criterion 
that gives the characteristic mass Mc at which the baryon 
fraction is on average half that of the cosmic mean. This 
is phrased on a condition on the virial temperature which 
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Figure 15. Stellar mass function at 2 = 0. We assign a stel- 
lar mass to eacli lialo in the following manner: if Mtot (2: = 0) > 
Mc(0) the stellar mass is (/*}f2(,/^m A^tot- For smaller halos we 
compute the redshift 2* at which Mtot exp(— 02*) = Mc{z*). We 
compute the stellar mass according to A/(2*). The 

solid line show the results for the halo sample from the basic run 
and the parameters (/*) = 0.8 and a = 0.35. Numbers along 
the line indicate how many halos are actually in our void sample. 
For comparison the stellar mass function obtained from the sim- 
ulation (dash-dotted) is shown. In addition, we show the result 
for the assumption that in each halo the stellar mass amounts to 
{f,)Ql,/^m Aftot (dashed line). 



must at least reach an "entry temperature" in the instable 
branch of the cooling/heating equilibrium line in the gas 
phase-space diagram. Our prediction for Mc (z) is strikingly 
well in accord with the evolution obtained from the simula- 
tions. 

At redshift 2 = 0, the characteristic mass scale of 
photo-evaporation corresponds to a circular velocity of ~ 
27 kms~^, which is significantly less than we wo uld ob- 
tain w ith the filtering mass formalism introduced bv lGnedlnl 
l|2000() using an entirely different approach. We thus pre- 
dict a considerably milder effect of photo-heating on the 
evolution of small halos. This has significant consequences 
for semi-analytic models of galaxy formation that relied on 
photo-heating as a feed back mechanism t o suppress small 
galaxies. For instance, ISomerviIl3 (|20o3) concluded that 
the substructure problem in the Local Group can be eas- 
ily solved taking photo-evaporation into account, provided 
the eff ect is as strong as suggested by the analysis of lGnedinl 
l|200(]|) . This conclusion is probably no longer valid when the 
effect is much weaker, as we find here. 

It is clear however that the baryon fraction and hence 
the star fraction rate in void dwarf galaxies is significantly 
reduced by the UV-background. As a result, the halo mass 
function cannot be translated directly into a mass function 
of the luminous matter. In the framework of the physics 
included in our simulations, we argue that the baryonic 
mass belonging to a halo is essentially set at the time when 
Mvii{z) — Mc{z). We compute the resulting mass function 



of the stellar content in halos under this assumption and 
obtain good agreement with our direct simulation measure- 
ments. 

We caution however that it is numerically difficult to 
resolve all the star formation in small dwarf galaxies at high 
redshift. While we have found some dwarf halos in our sim- 
ulation without any stars, it seems plausible that this is due 
to insufficient resolution, because all halos with a progenitor 
with M{z) > Mc{z) at some z should have been able to cre- 
ate some stars. Only in those halos with a very steep mass 
accretion history, steep enough to be for all redshifts below 
Mc{z), star formation should always be suppressed. 

Our results are derived under the assumption of op- 
tically thin gas with a spatially uniform UV fiux, neglect- 
ing any effects due to s clf-shiclding in the inner parts of 
halos. Recently. ISusa iST Umcmura (2004a) analysed photo- 
evaporation with 3D radiation transfer simulations for in- 
dividual halos. They found that for both a full radiation 
transfer treatment and a optical thin approximation the fast 
rising UV-background at high redshift efficiently suppresses 
star formation. However, a full radiation transfer treatment 
with self-shielding can only increase the amount of stars be- 
cause it makes the UV heating less efficient. Neglecting self- 
shielding in our analysis may thus lead to an overestimate 
of the effects of the UV-background. 

The evolution of the characteristic mass function de- 
pends essentially on the UV-fiux history. We have used 
here the Haardt & Madau ( 1996) model, where reionisation 
takes place at redsh ift z — 6. Previous results from WMAP 
JSoergel et al.ll2003l) suggested that the Universe could have 
been reionised at much higher redshift. If true, the evolution 
of Mc could be much more shallow at high redshift. Conse- 
quently, more halos would have an accretion history entirely 
below Mc{z), and hence more halos would remain virtually 
free of condensed baryons. The most recent analysis from 
three year data of WMAP iSoergel e t al. 2006) indicates 
that reionisation took place in the redshift rage z ~ 11 to 
7. This would lead to a moderately steeper Mc{z). Inter- 
estingly, the so modified characteristic mass function would 
evolve almost parallel to the average accretion history, al- 
lowing both, halos which condensate gas only at high red- 
shifts and others which start late to condensate. However, 
our measurements of the characteristic mass function based 
on the standard reionisation at 2 = 6 seem to be a good 
approximation for the effects of UV-heating in evaporating 
baryons of dwarf size halos. 

In summary, it appears unlikely that adjustments in 
the reionisation history can alter our basic finding that 
UV-heating is not particularly efficient in evaporating all 
baryons out of dwarf-sized halos. This makes it question- 
able whether feedback by photo-heating is really sufficient 
to suppress dwarf galaxies in voids, luminous satellites in 
galaxies, and to fiatten the faint end of the galaxy lumi- 
nosity function as much as observationally indicated. Other 
proposed solutions like kinetic supernova feedback and as- 
sociated galactic winds may therefore be needed to resolve 
these problems completely. 
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